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Abstract 

The paper deals with a lot sizing problem with ill-known demands modeled by fuzzy 
intervals whose membership functions are possibility distributions for the values of the 
uncertain demands. Optimization criteria, in the setting of possibility theory, that lead 
to choose robust production plans under fuzzy demands are given. Some algorithms 
for determining optimal robust production plans with respect to the proposed criteria, 
and for evaluating production plans are provided. Some computational experiments are 
presented. 
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1 Introduction 

Nowadays, companies do not compete as independent entities but as a part of collaborative 
supply chains. Uncertainty in demands creates a risk in a supply chain as backordering, 
obsolete inventory due to the bullwhip effect pQ . To reduce this risk two different approaches 
exist that are considered here. The first approach consists in a collaboration between the 
customer and the supplier and the second one consists in an integration of uncertainty into a 
planning process. 

The collaborative processes mainly aim to reduce a risk in a supply chain [2]. This is 
done by enforcing a coordination in a supply chain. Two approaches can be applied: vertical 
and horizontal. The vertical approach is a centralized decision making that synchronizes a 
supply chain (the most common way to coordinate within companies). The horizontal one 
refers to the collaborative planning, in which a supply chain can be seen as a chain, where ac- 
tors are independent entities [3]. The industrial collaborative planning has been standardized 
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for implementing a cooperation between retailers and manufactures. This process is called 
Collaborative Planning, Forecasting and Replenishment (CPFR) More precisely, the col- 
laborative processes are usually characterized by a set of point-to-point (customer/supplier) 
relationships with a partial information sharing [21 [5]. In the collaborative supply chain, a 
procurement plan is built and propagated through a supply chain. Namely, the procurement 
plan is composed of three horizons: freezing, flexible and free ones |2j. Quantities in the 
freezing horizon are crisp and can not be modified, quantities in the flexible horizon are in- 
tervals and can be modified under constraints imposed by a previous procurement plan. In 
the free horizon quantities can be modified without constraints. Another way to reduce a 
risk in a supply chain is to integrate the uncertainty in a planning process. In the literature, 
three different sources of uncertainty are distinguished (see j6] for a review): demand, process 
and supply. These uncertainties are due to difficulties to access to available historical data 
allowing to determine a probability distribution. 

In this paper, we focus on the collaborative supply chain (a supply chain, where actors 
are independents entities) under uncertain demands. In most companies today, especially 
in aeronautic companies, actors use the Manufacturing Resource Planning (MRPII) to plan 
theirs production. MRPII is a planning control process composed of three processes (the pro- 
duction process, the procurement process and the distribution process) and three levels [7J: 
the strategic level (computing commercial and industrial plans), the tactical level (the Master 
Production Scheduling (MPS) and the Material Requirement Planning (MRP)) and the oper- 
ational level (a detailed scheduling and a shop floor control). MRPII have been also extended 
to take into account: the imprecision on quantities of demands (MPS) [8J, the imprecision on 
quantities of demands and uncertain orders [9] (MRP) and the imprecision on quantities and 
on dates of demands with uncertain order dates [10J (MRP). 

In this paper, we wish to investigate the part of the MRPII process. Namely, the pro- 
curement process in the tactical level in the collaborative context. Our purpose is to help the 
decision maker of a procurement service to evaluate a performance of a given procurement 
plan with ill-known gross requirements and to compute a procurement plan in a collaborative 
supply chain (with and without supplier capacity sharing due to a procurement contract) 
with ill- known gross requirements. 

Several production planning problems have been adapted to the case of fuzzy demands: 
economic order quantity [TU [12], multi-period planning [8j [TUl El HH [151 [161 H7J> an d the 
problem of supply chain planning (production distribution, centralized supply chain) [18, 19, 
EQl EH [22l [23] . In the literature, there are two popular families of approaches for coping with 
fuzzy parameters. In the first family, a defuzzification is first performed and then deterministic 
optimization methods are used [20|, 12 lj . In the second one, the objective is expressed in the 
setting of possibility theory [21] and credibility theory [25]. We can distinguish: the possibilistic 
programming (a fuzzy mathematical programming) in which a solution optimizing a criterion 
based on the possibility measure is built |16,|T7], the credibility measure based programming 
in which the credibility measure is used to guaranty a service level (chance constraints on the 
inventory level) [26j or the goal is to choose a solution that optimizes a criterion based on the 
credibility measure [13] and a decision support based on the propagation of the uncertainty 
to the inventory level and backordering level [8j El [10]. Here, we restrict our attention to 
uncertainty propagation in MRP (the tactical level) [8j[9l[T0] and we propose methods both for 
evaluating a procurement plan in terms of costs under uncertain demands and for computing 
a procurement plan which minimizes the impact of uncertainty on costs, since the approaches 
proposed in the literature are not able to do this. 
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Popular setting of problems for hedging against uncertainty of parameters is robust op- 
timization |27| . In the robust optimization setting the uncertainty is modeled by specifying 
a set of all possible realizations of the parameters called scenarios. No probability distribu- 
tion in the scenario set is given. The value of each parameter may fall within a given closed 
interval and the set of scenarios is the Cartesian product of these intervals. Then, in order 
to choose a solution, two optimization criteria, called the min-max and the min-max regret, 
can be adopted. Under the min-max criterion, we seek a solution that minimizes the largest 
cost over all scenarios. Under the min-max regret criterion we wish to find a solution, which 
minimizes the largest deviation from optimum over all scenarios. 

In this paper, we are interested in computing a robust procurement plan (with and without 
delivering capacity of the supplier sharing). The delivering capacity are composed of two 
bounds: the lower one being the minimal accepted quantity that is sent to the customer 
and the upper bound which is due to a production capacity of the supplier. Moreover the 
customer accepts to have backordering but it is more penalized that inventory. This problem 
is equivalent to the problem of production planning with backordering, more precisely to a 
certain version of the lot sizing problem (see, e.g., [28, 29J), where: the procured quantities are 
production quantities, a production plan; delivering constraints are production constraints, 
capacity limits on production plans; and the gross requirements are demands. Thus, the 
problem consists in finding a production plan that fulfills capacity limits and minimizes the 
total cost of storage and backordering subject to the conditions of satisfying each demand. It 
is efficiently solvable when the demands are precisely known (see, e.g., [301 ED G-3 ) • However, 
the demands are seldom precisely known in advance and the uncertainty must be taken into 
account. 

In this paper, we consider the above problem with uncertain demands modeled by fuzzy 
intervals. The membership function of a fuzzy interval is a possibility distribution describing, 
for each value of the demand, the extent to which it is a possible value. In other words, 
it means that the value of this demand belongs to a A-cut of the fuzzy interval with the 
degree of necessity (confidence) 1 — A. To evaluate a production plan, we assign to it, degrees 
of possibility and necessity that its cost does not exceed a given threshold and a degree 
of necessity that costs of the plan fall within a given fuzzy goal. In order to find "robust 
solutions" under fuzzy demands, we apply two criteria. The first one consists in choosing a 
production plan which maximizes the degree of necessity (certainty) that its cost does not 
exceed a given threshold. The second criterion is weaker than the first one and consists in 
choosing a plan with the maximum degree of necessity that costs of the plan fall within a given 
fuzzy goal. A similar criterion has been proposed in [33] for discrete optimization problems 
with fuzzy costs. We provide some methods for finding a robust production plan with respect 
to the proposed criteria as well as for evaluating a given production plan under fuzzy-valued 
demands which heavy rely on methods for finding a robust production plan, called optimal 
robust production plan, in the problem of production planning under interval- valued demands 
with the robust min-max criterion. Namely, it turns out that the considered fuzzy problems 
can be reduced to examining a family of the interval problems with the min-max criterion. 
Therefore, we generalize in this way the min-max criterion under the interval structure of 
uncertainty to the fuzzy case. 

The paper is organized as follows. In Section [2j we recall some notions of possibility 
theory. In Section [3j we present a lot-size problem with backorders and precise demands. In 
Section [4j we present our results. Namely, we investigate the interval case, that is the lot-size 
problem with backorders in which uncertain demands are specified as closed intervals. We 
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construct algorithms for finding an optimal robust production plan (a polynomial algorithm 
for the case without capacity limits and an iterative algorithm for the case with capacity 
limits) and for evaluating a given production plan (linear and mixed integer programing 
methods, a pseudopolynomial algorithm). An experimental evidence of the efficiency of the 
proposed algorithms is provided. In Section[5j we extend our results from the previous section 
to the fuzzy case. We study the lot-size problem with backorders with uncertain demands 
modeled by fuzzy intervals in a setting of possibility theory. We provide methods for seeking a 
robust production plan with respect to two proposed criteria as well as for evaluating a given 
production plan under fuzzy- valued demands (the methods heavily rely on the ones from the 
interval case). The efficiency of the methods is confirmed experimentally. 

2 Selected Notions of Possibility Theory 

A fuzzy interval A is a fuzzy set in R whose membership function fj,^ is normal, quasi concave 
and upper semicontinuous. Usually, it is assumed that the support of a fuzzy interval is 
bounded. The main property of a fuzzy interval is the fact that all its A-cuts, that is the 
sets AW = {x : fi A {x) > A}, A G (0,1], are closed intervals. We will assume that A^ is the 
smallest closed set containing the support of A. So, every fuzzy interval A can be represented 
as a family of closed intervals A^ = [a~^, a + ^], parametrized by the value of A G [0, 1]. In 
many applications, the class of trapezoidal fuzzy intervals is used. A trapezoidal fuzzy interval, 
denoted by a quadruple A = (a, b, c, d) and its membership function has the following form: 

if z < a, 
if a < z < b, 

1 if b < z < c, 
J^c tfc<z<d, 
if z > d. 

Its A-cuts are simply [a + X(b — a),d — X(d — c)] for A G [0, 1]. Notice that this representation 
contains triangular fuzzy intervals (6 = c). 

Let us now recall the possibilistic interpretation of fuzzy intervals. Possibility theory |24j 
is an approach to handle incomplete information and it relies on two dual measures: possi- 
bility and necessity, which express plausibility and certainty of events. Both measures are 
built from a possibility distribution. Let a fuzzy interval A be attached with a single-valued 
variable a (uncertain real quantity). The membership function /j,-? is understood as a pos- 
sibility distribution, ir a = fj,-r, which describes the set of more or less plausible, mutually 
exclusive values of the variable a. It can encode a family of probability functions [33]. In 
particular, a degree of possibility can be viewed as the upper bound of a degree of proba- 
bility |34| . The value of n a (v) represents the possibility degree of the assignment a = v, i.e. 
Il(a = v) = TT a (v) = ^j(v), where Il(a = v) is the possibility of the event that a will take 
the value of v. In particular, ir a (v) = means that a = v is impossible and TT a (v) > means 
that a = v is plausible. Equivalently, it means that the value of a belongs to a A-cut A^ 
with confidence (or degree of necessity) 1 — A. A detailed interpretation of the possibility 
distribution and some methods of obtaining it from the possessed knowledge are described 
in [24} I35j . Let G be a fuzzy interval. Then "a G G" is a fuzzy event. The possibility of 
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"o G G", denoted by U(a G G), is as follows [36 : 

II(a G (?) = supmin{7r a (t;), (1) 

II(a G G) evaluates the extent to which "a G G" is possibly true. The necessity of event 
"a G G", denoted by N(a G G), is as follows: 

N(a G G) = 1 - n(a G) = 1 - supmin{vr a (t;), 1 - ^g(u)} (2) 

t>6R 

= inf max{l - TT a (v), fix (v)}, 

where 1 — [i^ is the membership function of the complement of the fuzzy set G. N(o G G) 
evaluates the extent to which "a G G" is certainly true. Observe that if G is a classical set, 
then Il(a G G) = sup^gg- K a (v) and N(a G G) = 1 — sup^^ G 7V a (v). 

3 The Deterministic Problem 

We are given T periods. For period t, t = 1, . . . ,T, let dt be the demand in period t, dt > 
(here we assume that the demands are precise), xt the production amount in period t, Xt > 0, 
Zt, u t the production capacity limits on xt- Let X C be the set of feasible production 
amounts. Two cases are distinguished, the case with no capacity limits X = {(x±, . . . ,Xt) '■ 
xt > 0, t = 1, . . . , T} and the one io#/t capacity limits X = {(xi, . . . , xt) '■ h < %t < ut,t = 
1, . . . , T}. Set D t = X^i=i an d = ^- =1 Xj, D 4 and X 4 stand for the cumulative demand 
up to period t and the production level up to period t, respectively. Obviously, X^_i < X t 
and D(_i < Dt, t = 2, . . . ,T. The costs of carrying one unit of inventory from period t to 
period t + 1 is given by c\ > and the costs of backordering one unit from period t + 1 to 
period t is given by cf > 0. The nonnegative real function Lt(u,v) represents either the cost 
of storing inventory from period t to period t + 1 or the cost of backordering quantity from 
period t + 1 to period t, namely L t (X t , Dj) = c£(Xf — D t ) if Xt > D t ; c^(D t — Xt) otherwise. 
The function has the form Lt(Xt, Dt) = maxjct (Xf — Dt), c^(Df — Xt)}. 

Our production planning problem with the deterministic (precise) demands consists in 
finding a feasible production plan x = (xi, . . . ,xt), x G X, that minimizes the total cost of 
storage and backordering subject to the conditions of satisfying each demand, namely 

T 

minF(^)=min^Lt(Xt,Dt). (3) 

t=l 

Obviously, the problem ^ is a version of the classical dynamic lot-size problem with backo- 
rders (see, e.g., [2SIES])- Without loss of generality, we can assume that an initial inventory Iq 
and an initial backorder Bq are equal to zero. Otherwise, one can append period and assign 
xq = Iq and do = with zero inventory cost if Iq > or assign xo = and do = Bq with zero 
backorder cost if Bq > 0. 

In the case with no capacity limits problem ([3]) has a trivial optimal solution equal to 
(d\, . . . ,dr)- hi the case with capacity limits, (f3p can be formulated as the minimum cost 
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flow problem (see, e.g., [3 



T 



min ^(cf/i + cf Si) 



t=i 

s.t. B t -It = Y, t j=1 (d j -x j ), t = l,...,T, ( 4 ) 

k < xt <u t , t = l,...,T, 

B t ,I t >0, t = l,...,T. 

Problem Q can efficiently solved, for instance, by an algorithm presented in [32] that takes 
into account a special structure of the underling network. 

4 Robust Problem 

Assume that demands dt, t= 1, . . . , T, in problem Q, are only known to belong to intervals 
Dt = [d^~,df ], dj > 0. This means that we neither know the exact demands, nor can we 
set them precisely. We assume that the demands are unrelated to one another and there is 
no probability distribution in D t , t = 1,...,T. A vector S = (d\, . . . ,dx), dt G D t , that 
represents an assignment of demands dt to periods t, t = 1, . . . , T, is called a scenario. Thus 
every scenario expresses a realization of the demands. We denote by T the set of all the 
scenarios, i.e. F = [c£[~, df] x • • • x [d^, dj,]. Among the scenarios of T extreme scenarios can 
be distinguished, that is the ones, which belong to {d± , d^ } x • • • x {d^, dj,}, the set of extreme 
scenarios is denoted by r ex t- We denote by S + (resp. S~) the extreme scenario in which all 
the demands are set to their upper (resp. lower) bounds. The demand and the cumulative 
demand in period t under scenario S are denoted by d t (S) G D t and T) t (S), respectively, 
D t (5) = E*=i MS)- Clearly, for every S G T it holds D t _i(S) < D t (5), t = 2, . . . , T, and 
D t (5) € [D t (5-),Di(5+)]. The function Lt(X t ,D t (5)) = max{cf(X t - D t (5)), cf (D t (5) - 
X()}, represents either the cost of storing inventory from period t to period t + 1 or the cost of 
backordering quantity from period t+1 to period t under scenario S. Now F(x, S) denotes the 
total cost of a production plan x G X under scenario S, i.e. F(x,S) = Ylt=i Lt(X-t, Dt(S)). 

In order to choose a robust production plan, one of robust criteria, called the min-max 
can be adopted (see, e.g. [27]). In the min-max version of problem Q, we seek a feasible 
production plan with the minimum the worst total cost over all scenarios, that is 

T 

ROB : min^4(a:) = min max F(x, S) = min max N Lt(Xj, Df (S)). 

In other words, we wish to find among all production plans the one that minimizes the 
maximum production plan cost over all scenarios, that minimizes A(x), A(x) is the maximal 
cost of production plan x. An optimal solution x r to the problem ROB is called optimal 
robust production plan. 

Let x G X be a given production plan. A scenario S° G T that minimizes the total 
cost F(x,S) of the production plan x is called optimistic scenario. A scenario S w G T that 
maximizes the total cost F(x, S) of the production plan x is called the worst case scenario. 

4.1 Evaluating Production Plan 

In this section, we show how to evaluate a given production plan x* G X. We first consider the 
problem of determining the optimal interval, F x * = [f~, , f^*], containing possible values of 
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costs of the production plan x* which can be rigorously defined as the following optimization 
problems: 



/- =minF(x*,S), (5) 

f±=maxF(x*,S). (6) 
x ser 

It is easily seen that the problem of computing the optimal lower bound on costs of x* 
([5]) is equivalent to the one of determining an optimistic scenario S° for x* , namely f~„ = 
F(x* , S°) = mins g r F(x* , S). Similarly, the problem of computing the optimal upper bound 
on costs of x* ^ is equivalent to the problem of determining a worst case scenario S w for 
x*, i.e. /+ = A(x*) = F(x*,S w ) = max<j er F(x*,S). Thus 

F x * = [f~*,m = [F(x*,S°),F(x*,S w )}. (7) 

Using the optimal interval F x * of possible values of costs of production plan x* allows 
us to evaluate possibility and necessity that the cost of the plan does not exceed a given 
threshold under uncertain demands modeled by intervals. Hence, in order to assert possibility 
that the cost of the plan does not exceed a given threshold g, i.e. to assert whether there 
exits a scenario S G T for which F(x*, S) < g, it suffices to determine an optimistic scenario 
S°, the optimal lower bound f~* = F(x*,S°) and evaluate f~* < g. If the inequality holds 
then there exists a scenario; otherwise not. Similarly, evaluating necessity that the cost of 
the plan does not exceed a given threshold g, i.e. asserting whether F(x*,S) < g for every 
scenario S G T, we only need to determine worst case scenario S w , the optimal upper bound 
f£* = F(x*,S w ) and evaluate , f£* < <?• Thus, evaluating a production plan boils down to 
computing its optimistic and worst case scenarios. 

Let us consider the problem of computing an optimistic scenario for a given production 
plan a;* G X, that is the problem ([5]). Its minimum is attained for some 5 £ T, since F(x, S) 
is a continuous function on the bounded closed set T. Problem §5§ can be formulated as a 
linear programming problem: 



/".= min Tl=Mh + cfB t ) 

s.t. B t -I t = Y J t j=l {s j -x*), t = l,...,T, 

s t e[dt,df], t = l,...,T, 

B t ,I t >0, t = l,...,T. 



(8) 



If s°, B° and 1°, t = 1, . . . ,T, is an optimal solution to problem ([8]), then S° = (s°, . . . , Sj) is 
an optimistic scenario for x* (an optimistic realization of uncertain demands) and If is storing 
inventory amount from period t to period t + 1 and B° represents backordering amount from 
period t + 1 to period t under the optimistic scenario S°. The problem ^ can be reduced to 
the classical minimum cost flow problem and effectively solved by algorithms that take into 
a special structure of the underling network (see, e.g., [32]). Hence, and fact that c[,cf > 
it follows that for t = 1, . . . , T one of 1° and B° is zero. 

Let us study the problem of computing a worst case scenario for a given production plan 
x* G X, that is the problem Since F(x*,S) is a continuous function on the bounded 
closed set T, it attains maximum for some S G T. The problem ^ can be formulated as a 
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ft 



mixed integer programming problem (MIP): 

max Y.^iicth + cfBt) 
s.t. B t -I t = Y? j= i{sj 
s t G [df,df], 
B t Jt>0, 



*), t = l, 
t = l, 
t = l, 

it<(i-St)Y,U( x *j- d J)> * = 1 > 
s t <^Ei=i(4-^)» * = i, 

^G{0,1}, i = l, 



T 
T 
T 
T 

,T, 
,T. 



(9) 



Let sf, Bf, If and <5(, t = 1, . . . ,T, be an optimal solution to problem ([9]). Then S w = 
(sf , . . . , Sy) is a worst case scenario for x* (a pessimistic realization of uncertain demands) 
and If is storing inventory amount from period t to period t + 1 and Bf is backordering 
amount from period i+1 to period i under the worst case scenario S" 1 ". The last two constraints 
model ([9]) and the binary variables 5t ensure that storing inventory from period t to period t+1 
and backordering from period t + 1 to period t is not performed simultaneously (either It > 
or i?( > 0). If #t = 1 then backordering is performed Bt > 0; otherwise storing inventory is 
performed It > 0. Thus, the problem ^ turns out to be much harder than ([5]). 

We now solve the problem ^ by means of dynamic programming. Let us present a result 
which shows that determining a worst case scenario S w can be restricted to the vertices of T, 
that is to the set of extreme scenarios r ex t- We prove the convexity of the cost function on T. 

Proposition 1. Function F(x*, S) is convex on T for any fixed production plan x* G X. 

Proof. Function c\ (XJ-D t (5)) and cf (D t (5)-XJ) are convex on T. Then so are max{c[ (X£ — 
D t (5U cf (D t (5) - X|)} and £f =1 max{c( (X* - D t (5)), cf (D t (5) - X*)}. □ 

The following result allows us to reduce the set of scenarios T to the set of extreme 
scenarios T ext . 

Proposition 2. An optimal scenario for problem |6p (a worst case scenario) is an extreme 
one. 

Proof. Function F(x*, S) attains its maximum in T. Since F(x*, S) is convex (Proposition [l]) 
and r is the hyper-rectangle, an optimal scenario for problem ^ is an extreme one (see, 
e.g., 138). □ 



Applying Proposition [2j we can rewrite problem ^ as: 
f+=A(x*) = F(x*,S w 



max F(x*,S). 



(10) 



We are now ready to give a dynamic programming based algorithm for solving problem ( 10 ). 
Let B>t be the set of feasible cumulative demand levels in period t, t = 1, . . . , T, i.e. B>t = 
{D t (S' _ ), T>t(S~) + 1, ... ,D t (S' + )}, let £(_i(D) be the maximal cost of a given production 
plan x* over periods t,...,T, when the cumulative demand level up to period £ — 1 is equal 
to D, D G D t _i, A-l : ^-1 -> K+- Set D = {0}. It is evident that: 



£t(D) = 
£t_i(D) = max 



L t (X?,D + ^)+A(D + d t -) 
L t (X*,D + d+)+A(D + d+) 



D G Or, (11) 
D G D t _i, (12) 
t = T, . . . , 1. 
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The maximal cost of production plan x* over period 1, . . . , T is equal to £q(0), A)(0) = / : 



X 



which is computed according to the backward recursion (11) and (12). The corresponding 



to x* worst case scenario S w can be determined by a forward recursion technique. It is 



sufficient to store for each D 6 the value for which the maximum in (12) is attained, 
that is either D + or D + df. The running time of the dynamic programming based 
algorithm is 0{T ■ D-r), which is pseudo-polynomial. We have thus proved the following 
theorem. 

Theorem 1. There is an algorithm for computing the maximal cost of a given production 
plan x* and its a worst case scenario S w , which runs in 0(T ■ Dj>). 

Since finding a worst case scenario requires taking into account only extreme demand 
scenarios (see Proposition [2]) , the running time of the above algorithm may be additionally 



refined by reducing the cardinalities of sets Of, t = 1, . . . , T, in (11 ) and (12). Note that we 
need only consider cumulative demand levels D which can be obtained by summing instanti- 
ated demands, at their lower or upper bounds, in the periods up to t. Namely, D G O t if and 
only if D = Y^\=i ^k, where dk € {d^,d^}. Hence, each reduced set of possible cumulative 
demand levels H>t has form {Dj, . . . , Dj}. Of course, I is bound by D^S^). An idea of the 
improved dynamic programming based algorithm can be outlined as follows: 
Step 1 (reducing sets JS) t ): built a directed acyclic network composed of node that is associated 
with Do = {0} and T layers that are associated with the reduced sets of possible cumulative 
demand levels B>t = {D|, . . . , D'}, t = 1, . . . ,T. The t-th layer is composed of nodes corre- 
sponding to cumulative demand levels D, D € Dj. The arc between D"_ x and exists if 
and only if either = D"_ x + d^ or = D"_ x + . Observe that, except for nodes from 
the last layer, each node in the network has exactly two outgoing arcs. An example of the 
constructed network is presented in Figure [T] 

Step 2: compute the maximal cost of production plan x* over period 1 , . . . , T according to the 



backward recursion (11) and (12) in the constructed network with the layers corresponding 



to reduced sets Dt and store for each D G Dt-i the value D £ D( for which the maximum in 



(12) is attained. 

Step 3: determine a worst case scenario for production plan x* by performing a simple forward 
recursion in the constructed network using the stored (in Step 2) values for which the maxima 



in (12) are attained. 

The network in Step 1 can be built in 0(T ■ max^i^,^ \^t\) time. The running time of 
Step 2 is the same time as Step 1. Step 3 can be done in 0(T) time. Hence, the overall running 
time of the improved algorithm is 0(T -max t= i y \^t\)- It is easily seen that maxj =1 x |Dt| 
is upper bounded by T)t(S + ) and at the worst case max t= i r ..^ |D 4 | = Dt(<S + ). 

Furthermore, the running time can be reduced if d'l — d^ = ■ ■ ■ = d£ — dTp = h. Then 
we find B t = {D t (S~), D t (5") + h, . . . ,T> t (S~) + th}, t = l,...,T. Now the running time is 
0(T 2 ), which is polynomial. 



4.2 Solving the Robust Problem 

Let us consider the problem ROB with no capacity limits, i.e. the problem with the set 
X = {(xi, . . . , xt) '■ Xt > 0, t = 1, . . . , T}. In this case, we make the assumption: the costs of 
carrying one unit of inventory from period t to period t + 1 for every t = 1, . . . , T are equal, 
we denote it by c 1 and the costs of backordering one unit from period t + 1 to period t for 
every t = 1, . . . , T are equal, we denote it by c B . Note that function F(x, S) is continuous on 
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Figure 1: An example of the constructed network in Step 1. 



X and T, T is a closed bounded set, and so A(x) is well defined continuous function on X (see, 
e.g., [391 Theorem 1.4]). We show that an optimal robust production plan x = (x±, . . . ,xt), 
i.e. x = arg min^x A(x), exists and can be computed by the following formulae: 



Xt == cBDt( ^ Dt(S ' ) , -n: X i X, .. t>2. 



(13) 



An algorithm for determining a production plan x according to ([13]) can be implemented in 
0(T) time. Before we show that x is optimal to problem ROB (an optimal robust production 
plan) with no capacity limits we prove the following proposition. 



Proposition 3. Let x be a plan computed according to (IS). Then x is feasible and S and 



S + are the worst case scenarios for plan x, i.e. A(x) = F(x, S~) = F(x, 5 + ). 
Proof. See Appendix [A] □ 
We are now ready to prove that x is an optimal robust production plan. 



Theorem 2. A production plan determined by formulae (13) is an optimal one for prob- 
lem ROB with no capacity limits. 

Proof. See Appendix [A] □ 

Note that if an initial backorder Bq or an initial inventory Jo are not equal to zero then 
one can modify the interval demand D\ as follows: D\ := [d^ + Bo,d^ + Bq] if Bq > or 



D\ := [d( — Iq, d±— Iq] if < Iq < dj , and apply formulae ( 13 ) to determine an optimal robust 
production plan for problem ROB with the modified demand. If Iq > d( then one appends 
period 0, as it has been described in Section [3j and applies an algorithm (Algorithm [TJ for 
the case with capacity limits, It = 0, 114 = M, t = 1, . .. ,T, where M is a large number. 

Let us turn to the problem ROB with capacity limits, i.e. the problem with the set 
X = {(xi, . . . , xt) ■ h < xt < 114, t = 1, . . . , T}. Notice T is a bounded closed set. Function 
F(x, S) is continuous on X and T and hence A(x) is continuous function on X (see, e.g., 
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Theorem 1.4]). From this and the fact X is a bounded closed set it follows that A(x) attains 
its minimum on X. 

We now construct an iterative algorithm for solving problem ROB based on on iterative 
relaxation scheme for min-max problems proposed in |40j . Similar methods were developed 
for min-max regret linear programming problems with an interval objective function |41|, 142] . 
Let us consider the problem (RX-ROB) being a relaxation of problem ROB that consists in 
replacing a given scenario set T with a discrete scenario set = {S 1 , . . . , S }, F^is C T: 

RX-ROB: a*=mina 

s.t. a>F(x,S k ) VS k £T dis , (14) 

iex, 

where S k = (s k )f =1 . The constraint a > F(x,S k ), called scenario cut, is associated with 
exactly one scenario S k G T^is- Since T^s Q I\ the maximal cost a* of an optimal solution 
x* of problem RX-ROB over discrete scenario set T^iis is a lower bound on the maximal cost 
of an optimal robust production plan x r for problem ROB, i.e. a* < A(x r ). Note that the 
scenario cut, a > F(x,S k ), associated with S k is not a linear constraint. One can linearize 
the cut by replacing it in RX-ROB with the following T + 1 constraints and 2T new decision 
variables: 

J TS k i „B r>S k \ 



Bf -if = £i=i(*}-*i), t = l,...,T, 
B? ,lf >0, t = l,...,T. 

Our algorithm (Algorithm [I]) starts with zero lower bound on the maximal cost of an 
optimal robust production plan x r , LB = 0, a candidate x* G X for an optimal solution for 
ROB and empty discrete scenario set, Tjis = 0- At each iteration, a worst case scenario 
S w for x* is computed by applying the method Q or the dynamic programming based 



algorithm presented in Section 4.1 Clearly, A(x*) = F(x*,S w ) is an upper bound on A(x r 
A(x r ) < A(x*). If a termination criterion is fulfilled (usually (A(x*) - LB)/ LB < e if LB > 1; 
A(x*) — LB < e otherwise, e > is a given tolerance) then algorithm stops with production 
plan x* , which is an approximation of an optimal robust production plan. Otherwise the 
worst case scenario S w = (sf)f =1 is added to I^is, the corresponding to S w scenario cut is 
appended to problem RX-ROB. Next the updated linear programming problem RX-ROB is 
solved to obtain a better candidate x* for an optimal solution for ROB and new lower bound 
LB = a*. Since set T^s is updated during the course of the algorithm, the computed values 
of lower bounds are nondecreasing sequence of their values. Then new iteration is started. 

In order to choose a good initial production plan x* G X in Algorithm [TJ we suggest 
to solve the classical production planning problem ^ with capacity limits (the model Q) 
under the midpoint demand scenario S m , i.e. dt(S m ) = {<L~ + df)/2, t = 1, . . . , T and take 
an optimal production plan under the midpoint scenario as an initial production plan. 

Theorem 3. Algorithm^ terminates in a finite number of steps for any given e > 0. 

Proof See Appendix [A] □ 

Note that if an initial inventory Iq or an initial backorder Bo are not equal to zero then 
one appends period 0, as it has been described in Section [3] and applies Algorithm [T] for T + 1 
periods. 
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Algorithm 1: Solving problem ROB. 



Input: Interval demands Dt = [d^~ , df], costs c{ , , t = 1, . . . T, initial production 

plan x* E X, a convergence tolerance parameter e > 0. 
Output: A production plan x r , an approximation of an optimal robust production 

plan, and its worst case scenario S w . 
Step 0. k : = 0, LB := 0, T dis : = 0. 
Step 1. x k := x*. 

Step 2. Compute a worst case scenario S w for x fc by applying the method ^ or the 
dynamic programming based algorithm presented in Section [4. 1| 
Step 3. A := F{x k , S w ) - LB. If LB > 1 then A := A/LB. 
If A < e then output x k , S w and STOP. 
Step 4. k:=k + l. 

Step 5. S k := S w , T dis := r dis U {S k } and append scenario cut a > F(x, S k ) to 
problem RX-ROB. 

Step 6. Compute an optimal solution (a;*, a*) for RX-ROB, LB := a*, and go to 
Step 1. 



Let us illustrate, by the following example, that solving problem ROB leads to a robust 
production plan. We are given 5 periods with the production capacity limits on a production 
plan: l x = 40, u\ = 50, l 2 = 30, u 2 = 40, l 3 = 30, u 3 = 40, l 4 = 10, u 4 = 35 and l 5 = 10, 
u§ = 35. The costs of carrying one unit of inventory from period t to period t + 1, cj-, 
for every t = 1 , . . . , 5 equal 1 and the costs of backordering one unit from period t + 1 to 
period t, cf , for every t = 1, ... ,5 equal 5. The knowledge about demands in each period 
is represented by the intervals: D x = [30,45], D 2 = [5,15], D 3 = [10,30], D 4 = [20,40] and 
D 5 = [20,40]. The scenario set F (states of the world) is V = [30,45] x [5,15] x [10,30] x 
[20,40] x [20,40] (see Figure [2]). The execution of Algorithm [l] gives a production plan: 
x° pt = 40, x° 2 pt = 30, x° 3 pt = 30, x° 4 pt = 27.9167, x° 5 pt = 10 with the total cost of 215.833 
(x opt is an approximation of an optimal robust production plan with convergence tolerance 
parameter e = 0.0001, the maximal cost of x op is no more than 0.01% from optimality). 
The worst case scenario S w G T is d^S™) = 30, d 2 (S w ) = 5, d 3 (S w ) = 10, d 4 (S w ) = 20, 
d 5 (S w ) = 20 and max Se r F{x opt , S) = F(x opt , S w ) = 215.833 (see Figure§). This means that 
total costs of production plan x opt do not exceed the value of 215.833 over the set of scenarios. 
Moreover, the plan x opt has the best worst performance, i.e. it minimizes the total cost over 



the all scenarios. Additionally, making use of the methods presented in Section 4.1, one gets 



complete information about all possible values of costs of the production plan x opt over the 
set of scenarios T, by determining the optimal interval F x o P t = [f~ opt , f^ op t] that contains 
these values (see Q). This interval equals [40,215.833]. A popular approach for solving a 
problem with uncertain parameters modeled by the classical intervals is taking the midpoints 
of the intervals (the average values of the possible parameter values) and solving the problem 
with these deterministic parameters. In our example, the midpoint scenario has the form 
dx(5 mid ) = 37.5, d 2 (S mid ) = 10, d 3 (S mid ) = 20, d 4 (S mid ) = 30, d 5 (S mid ) = 30. An algorithm 
for the problem with the midpoint demands (see Section [3]) returns an optimal production 
plan: x™ id = 40, xf id = 30, xf l = 30, x™ id = 10, x!f id = 17.5 with the total cost of 70. But, 
if scenario d^S™) = 45, d 2 (S w ) = 15, d 3 (S w ) = 30, d 4 (S w ) = 40, d 5 (S w ) = 40 (a worst case 
scenario for x mid ) occurs then the cost will be equal to 357.5 (max5 S r F(x mid , S) = 357.5). 
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Figure 2: The computed production plan x opt in the illustrative example. 



Table 1: Summary of the data and results of the illustrative example 
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robust 




plan 


plan 
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limits 


demands 


plan 


under S mid 


under S + 


under S~ 


t 


k u t 


A 


x° pt d t (S w ) 


x™ id d t (S w ) 


x+ d t {S w ) 


xt d t (S w ) 


1 


40 50 


[30,45] 


40 30 


40 


45 


45 30 


40 45 


2 


30 40 


[5,15] 


30 5 


30 


15 


30 5 


30 15 


3 


30 40 


[10,30] 


30 10 


30 


30 


30 10 


30 30 


4 


10 35 


[20,40] 


27.9167 20 


10 


40 


30 20 


10 40 


5 


10 35 


[20,40] 


10 20 


17.1 


40 


35 20 


10 40 


h 




= 5,t = 1, ...,5 


the worst costs for plans: 


„opt „mid _+ 








F(x opt ,S w ) 


F{x mia ,S w ) 


F{x+,S W ) 


F(x-,S W ) 








215.833 


357.5 


270 


395 



Note that 357.5 S> 215.833. Similar situation is for two extreme scenarios S + and S~ , the 
scenarios in which all the demands are set to their upper bounds and the lower bounds, 
respectively. Again running an algorithm for the crisp problem with the demands under 



scenario S + and S , we obtain optimal solutions: x{ 
xt = 35 with the total cost of 35 under S + and x7 



-- 45, x+ = 30, xj = 30, xj = 30, 
40, X2 = 30, X3 



30, x~l = 10, 



x 5 = 10, under S , with the total cost of 180. It turns out that if scenario ^(S"" 



30, 



d 2 (S w ) = 5, d 3 (S w ) = 10, d A (S w ) = 20, d 5 (S w ) = 20, S w G T, (a worst case scenario for x+) 
occurs then the cost of x + will be equal to 270. Similarly, if di(S w ) = 45, d 2 (S w ) = 15, 
d 3 (S w ) = 30, di(S w ) = 40, ^(S™) = 40, S w G F, (a worst case scenario for x~) occurs then 
the cost of x~ will be equal to 395. The summary of the input and output data of the above 
illustrative example is shown in Table [TJ Accordingly, we have no doubts that the computed 
plan x opt with respect to the min-max criterion (problem ROB) is a robust one. 

In order to check the efficiency of Algorithm [TJ we performed some computational tests. 
For every T = 100, 200, . . . , 1000, ten instances of the problem ROB with capacity lim- 
its were generated. In every instance, inventory costs were randomly chosen from the set 
{1,2,..., 10}, backorder costs were randomly chosen from the set {20, 21, ... , 50}, the de- 
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Table 2: Minimal, average and maximal computation times in seconds 



T 


100 


200 


300 


400 


500 


600 


700 


800 


900 


1000 


min 

avg 
max 


0.25 
0.39 
0.65 


0.63 
0.84 
1.05 


0.96 
5.68 
13.40 


2.04 
9.18 
20.47 


2.64 
11.47 
21.78 


23.55 
61.44 
147.32 


9.64 
72.20 
188.36 


24.54 
180.78 
328.61 


61.26 
282.20 
573.26 


41.44 
372.58 
893.01 



mands and the production capacities were randomly generated intervals [X, Y] , where X is 
an integer- valued random variable uniformly distributed in {0, 1, . . . , 99} and Y is an integer- 
valued random variable uniformly distributed in {100, 101, . . . , 199}. To solve the generated 
instances, we used IBM ILOG CPLEX 12.2 library (parallel using up 2 threads) |43j and 
a computer equipped with Intel Core 2 Duo 2.5 GHz. In Table [2j minimal, average and 
maximal computation times in seconds, required to find approximations of optimal robust 
production plans with convergence tolerance parameter e = 0.0001, are presented. Thus, the 
maximal costs of the computed production plans are no more than 0.01% from optimality. 
All computations finished in a few iterations and about 98% of the total running time was 
spent on computing worst case scenarios by MIP model Q. We also implemented and ran 
the improved dynamic programming based algorithm for computing worst case scenarios pre- 
sented in Section 4.1, but it turned out that solving MIP model for determining worst case 
scenarios was much faster than computing them by the dynamic programming algorithm. As 
we can see from the obtained results, Algorithm [I] allows us to solve quite large problems 
having up to 1000 periods in reasonable time. 



5 Fuzzy Problem 

In this section, we apply a more elaborate approach to model uncertain demands. Namely, 
the uncertain demands, in problem Q3p , are modeled by fuzzy intervals Dt, t = 1,...,T. 
Here, a membership function of Dt is regarded as a possibility distribution for the values of 
the unknown demand dt (see Section [2]). The possibility degree of the assignment dt = s is 
H(dt = s) = TTd t (s) = fig (s). Let S = (st)f =1 be a scenario that represents a state of the 
world where dt = st, for t = 1, . . . ,T. It is assumed that the demands are unrelated one 
to each other. Hence, the possibility distributions associated with the demands induce the 
following possibility distribution over all scenarios in S G M. T (see 



tt(S) = H((di = si) A • • • A [drp = st)) = mm^U(dt = St) = min /^g t (sj). (15) 

The value of ir(S) stands for the possibility of the event that scenario S 6l T will occur. We 
have thus extended scenario set T given by the intervals (see Section [4]) to the fuzzy case and 
now T is a fuzzy set of scenarios with membership function Hf(S) = vr(5), S G M?. We see 
at once that the A-cuts of T for every A G (0, 1] fulfill the following equality: 

fM = {S : tt(S) > A} = [d~ [X] ,df X] ] x ••• x [d T [X \dp X] ], 



which is from (15) and the definition of A-cut. We also define T>°> = [d l , d 1 ] x • • • x 
[d T ,di, ]. Notice that T A , A G [0,1], is the classical scenario set containing all scenarios 
whose possibility of occurrence is not less than A. 
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5.1 Evaluating Production Plan 

In order to choose a reasonable production plan under fuzziness, we first show how to evaluate 
a given production plan x G X. Notice that a cost of production plan x is unknown quantity, 
denoted by f x , since demands are unknown and modeled by fuzzy intervals in the setting of 
possibility theory. Thus, the unknown cost f x falls within fuzzy interval F x , called fuzzy cost 
of plan x, whose membership function /ip is a possibility distribution for the values of fuzzy 
variable f x (unknown cost of x), ir^ = fi^ , defined as follows: 

Hp (v) = n(/« =v)= sup tv(S), veR. (16) 

{S: F(x,S)=v} 



Making use of (16), we can define degrees of possibility and necessity that a cost of a given 



plan x G X does not exceed a given threshold g: 

n(/ x <g) = supvTp (v) = sup vr(5), (17) 

v<g X {S: F(x,S)<g} 

N(f x <g) = l-IL(f x >g) = l- S wp-K F (v) = l- sup tt(S). (18) 

i»3 ^ {5: F(as,S)>ff} 

It is easily seen that TL(f x < g) = A means that there exists a scenario S such that 7r(5) = A 
in which the cost of plan x does not exceed threshold g, F(x, S) < g. N(f x < g) = 1 — A 
means that for all scenarios S such that vr(S') > A, costs of plan x under these scenarios do 
not exceed g. 



We now consider the problem of computing the degrees (17) and (18) of a given plan x. 
Write fJ. A ' = [/aT^, fx^]- Note that the interval F X X ^ (A-cut of the fuzzy cost) is the optimal 
interval of possible costs of x (see Q) in problem ROB under interval scenario set rW. Hence 
there exists a link between the interval and the fuzzy cases: 

n(/« <g) = sup{A G [0, 1] : /" [A] < g}, (19) 
N(/«<s) = l-inf{AG [0,1] :/+ [A] <g}. (20) 



Prom equations (19) and (20), we obtain methods for computing the degrees. So, in order to 
compute H(f x < g) (resp. N(/a; < g)) we need to find the largest (resp. smallest) value of 
A such that there exits a scenario S G for which F(x, S) < g (resp. for every scenario 
S G r' A ' inequality F(x, S) < g holds), which is equivalent to determine an optimistic scenario 
S° G f W (resp. a worst case scenario S w G r^) for x by solving (IsJ) (resp. Km) and evaluating 



F(x,S°) < g (resp. F(x,S w ) < g). Notice f~ [X] = F(x,S°) (resp. /+ [A] = F(x,S w )). Since 
f x ^ (resp. f x ^) is nondecreasing (resp. nonincreasing) function of A, we can apply a binary 
search technique on A G [0, 1]. 

The fuzzy cost of production plan x (the possibility distribution for costs of a;), F x can be 
determined approximately, if necessary, via the use of A-cuts. Namely, the optimal intervals of 
possible costs of F x ^ = [f x ^, f x ^] under are computed for suitably chosen A-cuts. Then 
fuzzy cost F x is reconstructed from their A-cuts. This approach makes sense since intervals 
-Fj. A ' are nested. 

5.2 Fuzzy Robust Problem 

We now propose two criteria of choosing a robust solution in the fuzzy- valued problem ^ . 
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We are given a threshold g, and we would like to find a production plan which maximizes 
the degree of certainty (necessity) that its cost does not exceed threshold g. Thus, we would 
like to solve the following problem: 



maxN(/ x < g). 



(21) 



There are no doubts that an optimal production plan computed according to (21) is a robust 



one, since with the highest degree of certainty costs of the plan over scenarios will not exceed 



threshold g. By (20), it is easy to check that problem (21) is equivalent to the following 



mathematical programming problem: 



min A 
S.t. f: 



. [A1 < 9, 



Ae [0,1], 



(22) 



If A* is the optimal objective and x* is an optimal solution for problem (22 ) then N(f x * < g) = 
1 — A* and x* is an optimal production plan for (21 ). If (22) is infeasible then N(f x < g) = 
for all x GX. 



We now present a more general criterion of choosing a robust production plan than (21). 



Namely, suppose that a decision maker knows her/his preferences about a cost of a production 
plan f x and expresses it by a fuzzy goal G, which is a fuzzy interval with a bounded support 
and a nonincreasing upper semicontinuous membership function /ig : M — > [0,1] such that 
Hq(v) = 1 for v G [0,g]. The value of /ig(/x) is the extent to which cost f x of x satisfies the 
decision maker. Now the requirement " f x < <?" is replaced softer one, i.e. li f x G G" . So, 
by u2ft and (16) the necessity that event " f x G G" holds can be expressed as follows: 



N(/ x G G) 



1 - n(/ a G) 

1 - supmin{7r /a », 1 - ^(v)} 



(23) 



■(■en 



= l-supmin{7r(5),l-)u?;(F(x,5))}. 

s 

Thus, if N(/a; G G) = 1 — A means that for all scenarios S such that vr(S') > A, the degree 
that costs of plan x fall within fuzzy goal G, is not less than 1 — A. Note that N(/a; G G) is 
more general and weaker than N(f x < g). If ^q{v) = for v > g then they are the same. 
Moreover, N(/ a < g) < N(f x G G). 

Let us give the second criterion of choosing a robust plan. We are given a fuzzy goal G, 
and we wish to find a production plan which maximizes the necessity degree that costs of the 
plan fall within fuzzy goal G. Thus we need to solve the following optimization problem, 



maxNf/, G G). 



(24) 



problem, which is from (20) and (23): 



We check at once that problem ( 24 ) is equivalent to the following mathematical programming 



min A 

s.t. /+ [A] <3 +[1 - A] , 
AG [0,1], 
X G X. 



(25) 
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If (x*,\*) is an optimal solution for problem (25), then ~N(f x * G G) = 1 — A*. If (25) is 
infeasible then N(/a, G G) = for all a; G X. 



An algorithm for solving problem (25) (resp. (|22|)) is based on the standard binary search 
technique in [0, 1] (the interval of possible values of A) which follows from the fact that 
is nonincreasing and (resp. g) is nondecreasing function of A. We call the algorithm 

the binary search based algorithm. To find an optimal (x*,A*), x* G X, A* G [0,1], with a 
given error tolerance ^ > 0, we seek at each iteration, for a fixed A, a plan x G X satisfying 
fx^ < (jr+[ 1_A l (resp. fx < <?)> which boils down to seeking an optimal robust production 
plan x r and its worst case scenario under scenario set i.e. to solving problem ROB (see 
Section Q. Note that /+ [A] = A(x) and also that /+ [A] < g+l*-*] (resp. /+ [A] < g) holds 
for some X G X if and only if it holds for an optimal robust production plan under T^. 
Thus, at each iteration, we use either formulae (13) for the case without capacity limits or 
Algorithm [T] for the case with capacity limits. If the length of determined interval of possible 
values of A is less than £, then an optimal robust production plan x r for a fixed A, (x r ,X), 
is an approximation of an optimal solution for problem ( |25| ) (resp. (|22[)) with precision £. 
The running time of the above algorithm is 0(I(T) log time, where £ > is a given 
error tolerance and /(ITI) is time required for finding an optimal robust production plan and 
its worst case scenario under interval scenario set (the running time of either (13) or 
Algorithm [TJ . 

We now show, by the following illustrative example, that determining a production pla n 
maximizing the necessity degree that costs of the plan fall within fuzzy goal G (problem (24)) 
is a robust one under uncertain demands modeled by fuzzy intervals. We are given 5 periods 
with the production capacity limits on a production plan: l\ = 40, u\ = 50, I2 = 30, ui = 40, 
I3 = 30, 113 = 40, I4 = 10, U4 = 35 and I5 = 10, 115 = 35 and the same the costs of carrying 
one unit of inventory from period t to period t + 1, c\ = 1, t = 1,...,5, and the same the 
costs of backordering one unit from period t + 1 to period t ,cf = 5, t = 1, . . . , 5. The demand 
uncertainty in each period is represented by the triangular fuzzy intervals: D\ = (30, 37.5, 45), 
D 2 = (5,10,15), D 3 = (10,20,30), D 4 = (20,30,40) and D 5 = (20,30,40), regarded as 
possibility distributions for the values of the unknown demands. The fuzzy set of scenarios 
has the membership function: fJ-f(S) = ir(S) = min^i,...^ jUg (st), S G M 5 . The fuzzy 

goal G is trapezoidal fuzzy interval G = (0,0, 195.83,215.42), where 195.833 is the maximal 
cost of an optimal robust production plan for the problem ROB without capacity limits 
(an ideal supplier) and under the supports (the interval demands) of the fuzzy demands, 
i.e. D*f\ t = 1,...5. Thus, a production plan with the cost less than 195.833 is totally 
accepted and with the cost greater than 215.42 is not at all accepted. The binary search 
based algorithm outputs a production plan (with £ = 0.01 and e = 0.0001 for Algorithm [T]) : 
x° pt = 40, = 30, Xg Pt = 30, x opt = 25.3776, x^ pt = 10, that maximizes the necessity 
degree that costs of the plan fall within fuzzy goal G, N(/ x o P t G G) = 1 — A = 0.883. 
This means that for all scenarios S whose possibility of occurrence is greater than 0.117, 
7r(5) > A = 0.117, the degree of necessity (the degree of certainty) that total costs of plan x opt 
fall within fuzzy goal G, is not less than 0.883. Furthermore, we are sure that the total 
costs of the plan do not exceed 196.986 (the total cost at A-cut equal to 0.117) for every 
scenario S such that tt(S) > 0.117. We now apply existing approaches to our example. We 
first consider methods based on a defuzzification which take into account only one scenario 
resulting from a defuzzification of fuzzy parameters (demands) - see, e.g., |20l I21j . Applying, 
for instance, the index of Yager [35], we get crisp demands: di(S Y ) = 37.5, d2{S Y ) = 10, 
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Table 3: Summary of the data and results of the illustrative example 
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(30,37.5,45) 


40 
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40 
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30 
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30 


30 


30 


4 
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10 


10 
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10 35 
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17.5 
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= 5,t = 1,...,5 


the worst costs under S such that tt(S) > 0.117 for 


G 


= (0,0,195.83,215.42) 


x° pt 


x Y 
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196.986 


313.262 


313.262 
















N(/ x o Pt e G) 


N(/ X v e G) 


n(/ b bz e G) 








0.883 


0.593 


0.593 



d 3 (S Y ) = 20, d 4 (S Y ) = 30, d 5 (S Y ) = 30. An algorithm for the crisp dynamic lot-size problem 
with these demands (see Section K^) returns an optimal production plan: x Y = 40, x Y = 30, 
xj = 30, xj = 10, x Y = 17.5 with the total cost of 70. However, the cost of x Y may 
be even 313.262 for scenarios S such that tt(S) > 0.117, and so 313.262 S> 196.986. The 
necessity degree that costs of x Y fall within G equals 0.593, N(/ a .Y eG) = 0.593, which gives 
N(/ :e y € G) < N(f x o P t € G). Let us examine a possibilistic programming (a mathematical 
programming with fuzzy parameters), where solution concepts are based on the Bellman- 
Zadeh approach [36] - see, e.g., [161 El- I n this way, the assertion of the form "f x £ G", 
where f x is a cost of production plan x, is treated as a fuzzy constraint and values of the 
membership function \1q stand for degrees of satisfaction of this constraint - the fuzzy goal. 
In other words, the assertion induces Tl(fx £ G). The joint possibility distribution generated 
by the fuzzy goal as well as constraints with fuzzy demands is the minimum of the possibility 
of satisfaction of the fuzzy goal and the possibility of feasibility of the constraints and thus 
an optimal production plan is a plan, denoted by a; BZ , that maximizes the possibility degree 
of satisfaction both goal and the constraints. A trivial verification shows that for production 
plan x Y with the total cost of 70, the possibility degree of satisfaction the goal as well as 
the constraints is equal to 1 in our fuzzy problem under consideration - in problem ^ with 
triangular fuzzy demands Dt, t = 1, . . . , 5, and so a; BZ = x Y . However, as we have seen above, 
this plan is not a robust one. The summary of the data and results of the example is given 
in Table [H 

In order to evaluate the efficiency of the binary search based algorithm, we show some re- 
sults of computational experiments. For every number of periods T = 100, 200, . . . , 1000, ten 
instances of the problem (24) with capacity limits were generated. In every instance, inventory 
costs were randomly chosen from the set {1,2,..., 10}, backorder costs were randomly chosen 
from the set {20, 21, . . . , 50}, the production capacities were randomly generated intervals 
[X, Y], where X is an integer- valued random variable uniformly distributed in {0, 1, . . . , 99} 
and Y is an integer-valued random variable uniformly distributed in {100, 101, . . . , 199}, the 
demands are triangular fuzzy intervals with the supports equal to [0, 199] and the modal 
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Table 4: Average computation times in seconds 





/S 


T 


0.00 • c 


0.25 • c 


0.50 -c 


0.75 • c 


1.00 • c 


100 


0.96 


0.92 


0.88 


0.83 


0.92 


200 


2.36 


2.29 


2.26 


2.11 


2.10 


300 


7.22 


6.93 


6.81 


6.76 


6.86 


400 


14.68 


14.67 


14.55 


14.57 


14.55 


500 


27.51 


27.06 


26.60 


26.69 


26.64 


600 


46.73 


46.85 


46.72 


46.38 


46.49 


700 


111.72 


109.70 


109.41 


109.71 


109.43 


800 


223.22 


217.24 


217.89 


217.65 


217.85 


900 


243.51 


243.28 


243.87 


244.75 


245.21 


1000 


417.78 


416.86 


416.14 


414.14 


415.80 



values equal to Z, where Z is an integer-valued random variable uniformly distributed in 
{0, 1, . . . , 199}, the fuzzy goal G was modeled as a trapezoidal fuzzy interval G = (0, 0, c, d), 
where c was chosen as the maximal cost of an optimal robust production plan for the prob- 
lem ROB without capacity limits and under the interval demands being the supports of 
the generated triangular fuzzy demands. The values of d were equal to c + f3, for f3 £ 
{0.00 • c,0.25 • c,0.50 • c,0.75 • c, 1.00 • c} and the error tolerance f = 0.01. We used IBM 
ILOG CPLEX 12.2 library (parallel using up 2 threads) [13] and a computer equipped with 
Intel Core 2 Duo 2.5 GHz to solve the generated instances. In Table [4] average computation 
times in seconds are presented. As we can see from the obtained results, the binary search 
based algorithm, which calls Algorithm [T] at each iteration with convergence tolerance param- 
eter e = 0.0001, can solve efficiently the problem (24), with capacity limits, having up to 1000 
periods. 



6 Conclusion 

In this paper, we have proposed methods to compute a robust procurement plan in the 
collaborative supply chain, where the customer uses a version of MRP with ill-known demands 
to plan a production. This problem is a certain version of the lot sizing problem with ill- 
known demands modeled by fuzzy intervals, whose membership functions are regarded as 
possibility distributions for the values of the unknown demands. We have introduced, in this 
setting, the degrees of possibility and necessity that the cost of a plan does not exceed a 
given threshold and a degree of necessity that costs of a plan fall within a given fuzzy goal, 
which allows us to evaluate a given production plan. Moreover, we have provided methods 
for computing these degrees. For finding robust production plans under fuzzy demands, 
we have proposed two criteria: the first one consists in choosing a production plan which 
maximizes the degree of necessity that its cost does not exceed a given threshold, the second 
criterion is softer than the first one and consists in choosing a plan with the maximum degree 
of necessity that costs of the plan fall within a given fuzzy goal. We have constructed the 
algorithms for determining optimal robust production plans with respect to the criteria and 
confirmed their efficiency experimentally. The criteria are a generalization, to the fuzzy case, 
of the known from literature the min-max criterion. Consequently, we have shown in the 
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paper that there exists a link between interval uncertainty with the min-max criterion and 
possibilistic uncertainty with the necessity based criteria. It turns out that the evaluation of 
a production plan and choosing a plan in the fuzzy-valued problem are not harder than in the 
interval- valued case. The difficulty of solving the fuzzy problems lies in the interval case, since 
it is reduced to solving a small number of interval problems. Therefore, we have discussed 
first the interval-valued case. In this case, we have considered the problem of determining 
the optimal interval of possible costs of a production plan, which allowed us to evaluate the 
plan. Determining the optimal bounds of the interval boils down to computing optimistic 
and worst case scenarios. We have proposed linear programing based method for computing 
an optimistic scenario and mixed integer programing and dynamic programming methods for 
computing a worst case scenario. We have also identified a polynomial solvable case. For 
computing an optimal robust production plan, we have provided a polynomial algorithm and 
iterative one for the cases: with no capacity limits and with capacity limits, respectively. Then 
we have extended the methods introduced for the interval- valued problem to the fuzzy- valued 
one. 

There is still an open question concerning the complexity status of computing a worst case 
scenario of a given production plan. The problem is pseudopolynomially solvable and poly- 
nomially solvable under certain assumptions and seems to be a core of most of the problems 
considered in the paper. These assumptions are nearly realistic and make possible extension 
of our approach to the case where a procurement plan is given for a family of product. In 
other words, when the sum of quantities procured has to respect supplier capacity constraints 
which are computed from a previous procurement plan. This problem is equivalent to the 
multi-item capacitated lot sizing problem. The fact that the complexity status is still open 
creates the possibility to find a polynomial algorithm and to extend our approach to the 
multi-item, multi-level capacitated lot sizing problem without the assumptions. So, it is an 
interesting topic of further research. 



A Appendix 

Proposition^ It is easily seen that X t G \Dt(S~), D t (5+)], c 1 {± t - D t (5~)) = c B (D t (5+) 
X t ), t = 1,...,T, and xi > 0. Since D t _i(5~) < D t (5~) and D t _i(fl'+) < D t (S 4 
t > 2, ^ shows that x t > 0, t > 2. Hence L t (X. t ,'D t (S~)) = c 1 (X t 
Li(X t ,D t (5+)) = c B {B t (S+) - X t ), t = 1,...,T, and so F(x,S~) 

any scenario, S G T. Therefore, F(x,S) = ^maxfc 1 ^ - D t (5)), c B (B t (S) - X t )} < 
Er = imax{c 7 (Xi-D t (5-)), C B (Di(5+)-X t )} = £t=i L t (± t , B t (S-)) = ^ =1 L t (X tl D t (5+)) 

□ 



T> t (S~)) and 
F(x,S+). Let S be 



F(x,S~) = F(x,S+). 



Theorem^ We show that A(x) > A(x) for every i£l Consider any x G X. Let us modify 
x in the following way: 



X, 



D t (S~) ifX t <D t (5-) ) 
X; ifD t (5-)<Xi<D t (5+),t 
[D t (5+) ifX;>D t (S+). 



Now X' t ' G [D t (5-),Dt(5 + )]. From the feasibility of x , it follows that X^ < X.' t , t > 2. 
Hence and Dt-i(S~) < D t (5~) and D t _i(5+) < D t (5+), t > 2, we obtain X^ < X^' , t > 2 
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and, in consequence x" € X. Furthermore, it is easy to check that A(x') > A(x"). By the 
definition of A(x ), we have A(x ) > max{F(o; , S~),F(x , 5' + )}. We now only need to show 
that max{F(x" ,S~),F(x" ,S + )} > F(x,S~) = F(x,S+). 



Let us focus on a production plan x (determined by formulae (13)). Notice that function 
F{x, S~) = Y%=x c 7 (X t - Dt(S-)) (resp. F(x, S+) = £ t T =1 c B (F> t (S + ) - X t )) are the sum of 
linear and increasing (resp. decreasing) functions with respect to Xt, c/(Xf — D^S 1- )) > 
(resp. c B (F> t (S + ) - X t ) > 0) for X t 6 (D t (5 _ ),D t (5 + )] j * = 1,...,T. Hence, for each 
£ = 1, . . . , T there exists the intersection point in interval [Dt(S~), D<(S' + )]. It is easy to check 
thatX t G [D t (5-),Di(5+)] and c 1 (± t -T> t (S-)) = c B (D t (5+)-X t ), t = 1, . . . , T. Therefore, 
the points Xt, t = 1,...,T, are the intersection ones and max{F(x ,S~),F(x , S + )} > 
F(x,S~) = F(x,S+). Proposition [3] shows that F(x,S~) = F(x,S + ) = A(x) and we thus 
get A{x) > A(x") > A(x). □ 

Theorem^ The proof is almost the same as those given in |47[ Theorem 2.5] and [4Q|, The- 
orem 3]. We will denote by {(x k ,a k )} the sequence of optimal solution (x*,a*) computed in 
consecutive iterations in Step 6 of Algorithm [TJ k stands for the k-th iteration. By picking 
a subsequence, if necessary, the sequence {(x k ,a k )} converges to point (x,a), X € X, which 
follows from the fact that sequence {x k } belong to the bounded and closet set X C (X is 
a compact set) and {a k } is a nondecreasing sequence bounded above. Similar considerations 
apply to the sequence {S k } of worst case scenarios determined in Step 2 of Algorithm [!} 
The set F is a closed and bounded (compact) and hence {S k } converges to S € F. Since 
scenario cuts are appended to problem RX-ROB, the inequality a k+l > F(x k+1 , S k ) holds. 
By continuity of F, we have 

a>F(x,S). (26) 

Let us define the set S w (x) of worst case scenarios for a; G X, i.e. S w (x) = {S w \ S w = 
arg maxggr F(x, S)}, which is the point-to-set mapping. The set S w (x) is nonempty for 
every for x £ X. By [391 Theorem 1.5], S w is upper semicontinuous at x and so S £ S w (x). 
Therefore 

A(x)=maxF(x,S) = F(x,§). (27) 



Combining (26) and (27), we obtain a > A(x). By [391 Lemma 1.2] A is upper semicontinuous, 
which yields 

A(x k ) = F(x k ,S k ) <a k + e, for some sufficiently large k. 

This implies that the termination criterion in Step 3 of Algorithm [T] will be fulfilled in a finite 
number of iterations. □ 
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